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ABSTRACT 

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local 
Universe (12 + log(0/II) — 7.17). In this work we study I Zw 18 using data from Spitzer, Herschel 
Space Telescope and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive 
maps of I Zw 18, to date, in both, the far infrared and the CO J = 1 transition. We use dust 
emission models to derive a dust mass upper limit of only Mdust < 1.1 x 10"* M© (3cr limit). This 
upper limit is driven by the non-detection at 160 /im, and it is a factor of 4 — 10 times smaller than 
previous estimates (depending upon the model used). We also estimate an upper limit to the total 
dust-to-gas mass ratio of Moust/^gas < 5.0 x 10~^. If a linear correlation between the dust-to-gas 
mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9x10^^. We 
also show that the infrared SED is similar to that of starbursting systems. 
Subject headings: galaxies: dwarf — galaxies: ISM — galaxies: individual(I Zw 18) 



1. INTRODUCTION 

The link between dust-to-gas mass ratio (DGR) and 
heavy element abundance ( metallicity) in galax i es re- 
mams an open issue (e.g iLisenfeld fc Ferraral 119981: 
[Ed mun 

dlEool iHunt et ahl I2005D . Specifically, in very 
low metallicity systems (12-|-log(0/II) < 8) it is unclear 
how the DGR scales with metallici ty. Models considering 
dust destruction by supernovae (iHirashita et al.l I2002D 
ormass outflows from the galaxv fcisenfeld fc Ferraral 
Il998f ) predict a nonlinear relation. On the other hand, if 
the fraction of meta ls incorporated in the dust is constant 
pames et al.ll200^ . we expect a linear relation between 
DGR and metallicity, in a sense that the ratio decreases 
as metallicity decreases. Measurements of DGRs over a 
range of metallicity are necessary to better constrain this 
relationship. 

The blue compact dwarf galaxy I Zw 18 has one 
of the lowest nebular met allicities measured to date. 
ISkillman fc Ke nnicutt (1993) measure an oxygen abun- 
dance of 12-hlog(0/H) = 7.1 7. This is 3.2% of th e solar 
abundance (using the scale of lAsplund et al.l[2009| ) . Most 
local universe galaxies have 12-|-log(0/II) ~ 8.5 (e.g., 
IMoustakas eFalllMol for S INGS), and the Milky Way 
has 12-Hog(0/II)~ 8.7 (jBaumgartner fc Mushotzkvl 
l2006f l. I Zw 18 therefore represents the extreme low end 
of the metallicity range in the local universe and is thus 
a key datum for understanding the relationship between 
DGR and metallicity. 

The dust mass of I Zw 18 is poorly known. Typi- 
cal galaxies of similar morphology (blue compact dwarfs) 
have dust masses that range between 10'^ — 10^ M©, with 

' Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA. 

^ National Radio Astronomy Observatory, 520 Edgemont 
Road, Charlottesville, VA 22903-2475, USA. 

^ Max-Planck-Institut fiir Astronomic, Konigstiihl 17, D- 
69117 Heidelberg, Germany. 

* Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218. 

^ Department of Physics & Astronomy, Macalester College, 
1600 Grand Avenue, Saint Paul, MN 55105, USA. 



DGR ranging between 10 — 10 ^ (jLisenfeld fc Ferraral 
19981). Usin g Ha/II/3 flux ratios as a dust tracer, 



Cannon et al.l (2002) find a total dust mass for I Zw 18 
of ~ 7 - 10 X 10^ Mel!| by assuming a linear scaling 
between DGR an d meta llicity (as measured by O/H). 
lEngelbracht et al.l (|2008f ). using Spitzer data limited by 
a non-detection at 160 /xm, measure an upper limit for 
the dust mass of 4.2 x 10^ Mp,. A more recent study of 
a large sample bv lGalametz et all ()2011[ ) uses previously 
published Spitzer and SCUBA data to constrain the dust 
mass of I Zw 18 to be < 1.1 x 10^ Mq and the DGR to 
< 4.5 X 10-4. 

I Zw 18 contains intense radiation fields stemming 
from active star formation. It therefore provides a 
nearby testing ground to probe the physics of dis- 
tant primeval sources. Previous studies using Spitzer 
(jEngelbracht et al.l 120081 : iWu et al.l 120071 ) show that its 
continuum emission from 15 to 70 iim has a slope char- 
acteristic of a starburst galaxy of solar abundance. More- 
over, its mid-infrared spectrum from 5 to 36 /im shows 
no detectable emission from polycyclic aromatic hydro- 
carbons (PAHs) . Such low abundance of PAHs is likely 
the consequence of a high radiation field in combination 
with the low metallicity of the source. 

In this paper we estimate the DGR for I Zw 18. We 
use previously unpublished Spitzer Space Telescope and 
archival Herschel Space Observatoryj continuum obser- 
vatio ns, combined \yith d ust emission models and a gas 
mass (|van Zee et al.|[l998l ) to constrain the radiation field 
intensity, temperature, dust mass and DGR in I Zw 18. 
Throughout this paper we assume a distance of 18.2 Mpc 
(Aloisi et al. 2007i ). Revisions to this distance will affect 



6 We scale ICannon et"all (200^ and 'Engclbracht et al.' 
result by a factor of (18. 2/12. 6)^ and|Galamctz et al. ( 2011) result 
by a factor of (18.2/13)^ to account for the differences in assumed 
distances. We note that our final result, the DGR, is distance 
independent. 

^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and v^ith 
important participation from NASA. 
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our dust and gas masses, but not the DGR. 

2. METHODS 
2.1. Observations 

We use a variety of data from several different fa- 
cilities to map the far infrared, sub millimeter, mil- 
limeter and radio wave emission of I Zw 18. There- 
fore our data set uses the following observatories and 
instrun ients: Spitzer M ultiband Imaging Spectrometer 
(MIPS. lRieke et al. 1'2004'): 77ersc/ieZ Photodctcctor Array 
Camera and Spectrometer (PACS, Poglitsch et al. 2010); 
Herschel Spectral and Ph otometric Imaging Receiver 
fSPIRE. [Griffin et al.llMol ): IRAM Plateau de Bure In- 
terferometer (PdBI) and Very Large Array (VLA). In 
this section we will briefly describe new observations. 

Spitzer: We observed I Zw 18 at 70 and 160 /im 
using MIPS in photometry mode as part of a cycle 5 
proposal (P.I. A. Bolatto, AOR: 22369536). The total 
observation time was 8 hours. The reduction of these 
images very closely follo ws the procedure descr ibed in 
iGordon et all (|2007| ) and lStansberrv et all ()2007f ) for the 
70 and 160 fim map respectively. 

Herschel /PACS: We use archival 70 and 160 /im ob- 
servations from Herschel . The observations were taken 
with PACS using the Large Scan Map mode as part of the 
Herschel Guaranteed Time Key Program, Dwarf Galaxy 
Survey (P.I. S. Madden, obs. ID: 1342187135/36). The 
scan maps were taken at 90° angles from one another at 
the medium scan speed (20 " s~^) and then combined 
together in order to reduce the noise caused by streaking 
along the scan direction. The scan leg length is 4.0 ' and 
the total on source time for the combined images was 
192 s. 

Unlike Spitzer data, the methods to reduce PACS data 
are still evolving significantly. Therefore, we reduce the 
data in two separate ways. We first use Herschel Interac- 
tive Processing Environment (HIPE) v4.2 with the stan- 
dard pipeline scripts. We also process the data up to level 
1 in HIPE v7. We use the standard pipeline, which in- 
cludes pixel flagging, flux density conversion, and sky co- 
ordinate association for each pixel of the detector. At this 
stage, the PACS timelines are still affected by 1// noise 
and baseline drifts. In order to subtract the baseline, 
remove glitches, and project the timelines on the final 
map, we applied the scanamorphos algorithm (jRoussell 
120121 ) to the levcl-1 PACS timelines. 

Herschel /SPIRE: We use archival Her schel Spectral 
and P hotometric Imaging Receiver (SPIRE: iGriffin et al.l 
120101) photometric observations at 250, 350, and 500 /xm. 
Observations were made in the large map mode with 
the nominal scan speed of 30 arcsec s~^ and the cross- 
scanning method as part of Herschel Science Demonstra- 
tion Phase (PI: S. Madden; obs. ID 1342188663). Data 
reprocessing was carried out in HIPE using the standard 
large map pipeline with the latest SPIRE calibration tree 
availablco which includes deglitching the timeline data, 
flux calibration and various corrections. After removal of 
a linear baseline, images were made using the standard 
naive mapper. The final maps are in units Jy beam"^ 
with pixel scales of 6, 10 and 14 arcsec at 250, 350 and 
500 fim, respectively as described in the SPIRE Data 
Reduction Guide. 

* We used HIPE 8.1 and the SPIRE calibration tree v. 8.1. 



PdBI: We present new observations of the CO J = 
1^0 transition in I Zw 18 using the IRAM Plateau 
de Bure Interferometer as project t027 (P.I. A. Leroy). 
The data were observed on 24, 27, and 28 September 
2009 using the " 5Dq" configuration, meaning that 5 tele- 
scopes were operational and that the array was in a 
compact configuration. They data is calibrated in the 
standard way in December 2009 using the PdBI pipeline 
implemented in the CLIC and MAPPING packages of 
GILDAS. The effective time on source was 12.5 hours af- 
ter fiagging during the pipeline run. The effective band- 
width was ^ 850 MHz, or about 2200 km s^^ with native 
resolution ~ 2.5 MHz (6.5 km s~^). We do not detect CO 
emission. At 26 km velocity resolution we achieved 
an RMS noise of 1.26 mJy beam~^, implying a 4cr flux 
upper limit for a point source of 0.131 Jy km s~^. 

VLA: The obse rvations used to cons truct the HI map 
are described in Ivan Zee et al.l (|1998[ ). We obtained 
two hours of Rapid Response 21 cm VLA observations 
(project 08B-246; P.I. A. Bolatto) to evaluate the Galac- 
tic foreground contribution. This contribution can be 
estimated by measuring the HI column density toward 
I Zw 18 and convert it to dust emission using t ypical 
high-latitude Galactic ratios (e.g. lBoulanger et al.i ri996'). 
The observations were obtained during the move between 
CnD and D configuration with a synthesized beam size of 
67" x41", and a native resolution of 6.1 kHz (1.3 km s~^). 
Data was reduced in AIPS using the standard procedure 
and calibrations, and care was taken to remove the base- 
lines affected by the frequency aliasing problems due to 
the VLA-JVLA transition. At 10.3 km s~^ velocity res- 
olution we achieved an RMS noise of 1.1 mJy beam~^, 
implying an HI column density of Nhi — 2.4 x 10^^ cm~^. 
Galactic neutral hydrogen emission was observed in the 
central 30 km s~^ of the passband. Nonetheless, even 
after spatially filtering the 160 /xm MIPS map to ap- 
proximately match the uv coverage of the VLA the cor- 
relation between the high resolution HI column density 
and the 160 /im surface brightness remained extremely 
low (Fig. [5]) , showing that most of the emission present 
in the 160 and 250 /im images is not due to the high- 
latitude Galactic foreground. 

2.2. Photometry 

In Fig. [l]we present the Spitzer MIPS 70 and 160 /im 
maps and the Herschel PACS 160 /im image for I Zw 18. 
Overlaid on each image is the HI column density distribu- 
tion observed by Ivan Zee et al.l fl998). The HI contours 
correspond to 0.7, 1.4 and 5 x 10^" cm~^ enclosing 98, 
96 and 78% of the total flux at 70 /im. We use the 
MIPS 70 /im map over the PACS 70 /im because it has 
a much better surface brightness sensitivity (0.17 versus 
2.87 MJy sr"-'^), yielding a better signal to noise. The 
bulk of the 70 /im emission coincides with the HI col- 
umn density maximum. This peak also coincides with 
the location o f activ e star forming regions observed by 
'Cann on et all (|2002h using the Hubble Space Telescope 
(HST). The diffuse component at 70 /im extends prefer- 
entially ~3 kpc north-west from the peak. We subtract 
the background emission measured in a region free of 
sources. We integrate the flux using a circular aperture 
with radius of 45" centered at the peak of the 70 /im emis- 
sion and applying an aperture correction factor of 1.13 
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Fig. 1. — The left two panels show the Spitzer images of I Zw 18 at 70 and 160 iim. The right pane l shows the Herschel PACS image 
at 160 fim. Overlaid as black contours is the Hi column density distribution from Ivan Zee et al.l ll 19981) using the VLA. The contours are 
0.7, 1.4 and 5 X lO'^" cm~^. The black circle in the bottom right corner of each panel corresponds to the respective beam size of the FIR 
observations. The smaller circle in the first panel corresponds to the beam size of the Hi observations. At 70 fj,m, the bulk of the emission 
coincides with the Hi contours and the diffuse emission extends preferentially towards the NW. I Zw 18 is undetected in both 160 fim 
maps. 
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Fig. 2. — Spitzer 160 fim (left) and Herschel 250 fim (right) maps of a ~5 X 5 arcmin field around I Zw 18. The thick black contours 
show the VLA observations of the Galactic Hi foreground emission at 2, 4 and 6cr significance level. The thin black contours are the same 
as shown in Fig. [T] The VLA beam size is 9 = 67.3" X 41.3" and is shown in the corner of the left panel. The SPIRE 250 ^m beam is 
shown in the right panel. The bulk of the 160 fim emission, located south-east of our object, fragments into at least three point sources 
in the 250 fim map, which has better spatial resolution. There is no correlation between the far infrared and the Galactic Hi foreground 
emission, suggesting that the confusion is dominated by background galaxies. 
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(determined by integrating over the point spread func- 
tion and compatible with those in the MIPS Instrument 
Handbook). The cahbration error on Spitzer is about 5% 
at 70/xm. We estimate a photometry error of 1.7 mJy 
by adding in quadrature the calibration uncertainty and 
the background noise We measure a total flux density of 
33.6±1.7 mJy at 70 /im. Ou r flux value is consistent w ith 
the 70 ^m flux measured by lEngelbracht et all ()2008[ ) of 
34.9±4.79 mJy. 

The Spitzer 160 fim map is confusion-limited. The 
bulk of the emission is associated with sources outside 
the HI emitting region of I Zw 18. Although the Galac- 
tic latitude of I Zw 18 is '^44°, it could be possible that a 
significant source of confusion were Galactic cirrus. We 
use the VLA HI observations to explore this possibility. 
Fig. El shows the MIPS 160 fim and the SPIRE 250 fim 
map of a ^5x5 arcmin field around I Zw 18. The thick 
lines represent the HI foreground emission from VLA. 
The thin lines represent the HI emission of I Zw 18. It is 
clear from visual inspection that the maxima of the HI 
foreground and the 160 /im and 250 /im emission are not 
coincident. We find a Pearson correlation coefficient close 
to zero. In most of the Spitzer confusion-limited images 
at 160 /im the confusion i s mainly due to faint unre- 
solved background sources (|Dole et al.1 12001 . The bulk 
of the 160 /im emission, located south-east of our object, 
fragments into at least three point sources in the 250 /im 
map, which has finer spatial resolution. The 160 /im peak 
also coincides with several background galaxies in deep B 
and i?-band images (S. Janowiecki, private communica- 
tion). The difference between the 160 and 250 /xm maps 
is consistent with what one would observe if the peak of 
the emission at a 160 /im is associated with background 
galaxies. 

The background contamination and the absence of cor- 
relation between the 160/250 /im emission and the HI 
foreground emission makes it impossible to recover the 
flux associated with I Zw 18. Thus, we use an annular 
sector around the galaxy to measure a one-sigma surface 
brightness sensitivity of 0.18 mJy sr~^ that includes the 
effects of confusion. We estimate the flux upper limit 
multiplying this value by the area associated with the 
1.4x10 atoms cm~^ HI contour that encloses 96% of 
the 70 /xm flux and 62% of the HI mass. To find the 
aperture correction factor associated to this area, we 
can approximate the contour using a circular aperture 
of 48" radius. This corresponds to an aperture correc- 
tion factor of 1.6 at 160 /im. After applying the aperture 
correction factor, we obtain a corresponding Str flux up- 
per limit of 40.5 mJy. This new upper limit is a factor 
of ^2 lower than t he pre vious upper limit published by 
lEngelbrachtltall (pOOSi ). 

The Herschel PACS image at 160 /im also fails to de- 
tect I Zw 18. In this case, however, the image is not 
confusion-limited. The la surface brightness sensitivity 
is 1.8 MJy sr~^. If we assume that the emission from 
I Zw 18 is compact on 12" scales, the corresponding 3cr 
flux upper limit integrating over the 12" beam and ap- 
plying an aperture correction factor of 1.32 is 27.2 mJy. 

We will work with the PACS 160 /im flux upper limit 
of 27.2 mJy for the rest of the paper. However, the up- 
per limit measured from the PAGS data relies upon the 
assumption that the source is compact. If the 160 /im 
emission is significantly extended over scales larger than 
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Fig. 3. — A Infrared SED of I Zw 18. Open symbols represent 3cr 
upper limits. The circles show the new MIPS data. The squares 
show the PACS and SPIRE data. The open diamond and trian- 
gle at 16 and 850 /im corresponds t o upper l imits estimated by 
lEngelbr acht et al. (2008) and Galame tz et al.l I I2011I) respectively. 
The solid line corresponds to the modified blackbody fit to the 
MIPS 70 fim detection and the PACS 160 fim upper limit. 

12" (~1 kpc), it may be more appropriate to use the 
MIPS 160 /im upper limit of 40.5 mJy. 

Finally, the SPIRE maps at 250, 350 and 500 /im show 
no detection of I Zw 18. From these images we measure 
a surface brightness sensitivity using an annular sector 
around the source. We then apply aperture corrections 
and point source color corrections assuming /3 = 1.5 (/3 in 
fi, (X i/^) described in the SPIRE Photometry Cookbook 
(Bendo, G. J. and the SPIRE-ICC, 2011) We measure 
3cr fiux upper limits of 22.2, 23.9 and 25.4 mJy at 250, 
350 and 500 /im, respectively. 

3. RESULTS 

3.1. Dust Mass 

We use two methods to estimate the dust mass of 
I Z w 18. In the first , we follow the procedure outlined 
bv iHildebrandl ()1983[ ). assuming an idealized graybody 
source with a single temp erature. In the second one we 
use the lDraine fcll ()2007D (DL07) model. The main dif- 
ference between the DL07 model and the idealized gray- 
body is that DL07 assumes a grain size distribution that 
reproduces the observed wavelength-dependence extinc- 
tion in the Milky Way, and consequently a distribution 
of temperatures. Giving the extreme nature of I Zw 18, 
it is not clear that either model is exactly applicable. 
Nonetheless, we used them so we can make a consistent 
comparison to larger samples of galaxies. 

3.1.1. Modified Blackbody Model 

For an idealized cloud, the dust mass is estimated 
by fitting its far infrared spectrum as the product of 
a blackbody spectrum {Bx,t) and a mass absorption 
coefficient (kx). The absorption coefficient varies with 
wavelength as the negative power of the grain emissivity 
index (ka oc A"'', where (3 represents the emissivity in- 
dex) . Then, for a cloud that is optically thick to starlight 
and optically thin to far infrared emission, the dust mass 
^Dust is given by the following expression: 

F 

Moust^^^ , (1) 
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where D is the distance to the galaxy, Fa is the observed 
flux and Sa.t is the blackbody intensity. The flux at 
any point on a graybody spectrum is Fa oc -Ba.tA"^, 
with /? independent of wavelength. Therefore, we can 
solve for the color temperature (T70/160) using the ratio 
F70 / F160. 

Measured value s for kx at 250 /xm (K250) span the range 
« 5 — 15 cm^ g~^ ()Alton et al.|[2004D . and commonly used 
values for (3 are 1 — 2 depending on the environment. For 
this work we adopt K250 = 9.5 cm^ g~^ and /3 = 1.5; 
these a re commonly used v alues for low metallicity galax- 
ies (e.g. |Lerox!iE3]llQQl3)- Using /3 = 1 or 2 changes our 
dust mass limits by ^ 10%. 

Fig. [3] shows the spectral energy distribution (SED) 
of I Zw 18. The source is only detected at 70 /im. At 
longer wavelengths each point corresponds to a Scr flux 
upper limits. Among these limits, the 160 /ini upper 
limit represents the strongest constraint on the I Zw 18 
SED. Thus, based on the 70 and 160 /im emission, the 
modified blackbody spectrum model constrains the dust 
temperature to be T70/160 > 33.7 K. This translates into 
a predicted 850 /im flux of 0.28 m Jy, compatible with th e 
observed upper limit of 1.25 mJy (jGalametz et al.|[2011[ ). 
Combining our temperature lower limit with the 70 /xm 
flux, Eq.[T]yields a dust mass oiMoust < 3.2 x 10^ M©. 
If we use the MIPS 160 fim upper limit instead of 
the PACS upper limit, we measure a temperature limit 
T70/160 > 29.8 K and a dust mass a factor of r^2 higher, 
i.e., Moust < 6.9 X 103 Mq 

3.1.2. Drame & Li Model 

For a det a iled d escription we refer to DL07 and 
IDraine et al.1 (|2007D . Essentially, DL07 models charac- 
terize the dust as a mixture of carbonaceous and amor- 
phous silicate grains with size distributions chosen to 
match the observed extinction in the Milky Way. To 
characterize the intensity of the radiation that is heat- 
ing the dust, the model adopts the spectrum of the local 
interstellar radiation field (this may not be a good ap- 
proximation for I Zw 18, a starburst system characterized 
by high intensity radiation fields and low metallicity). In 
DL07, most dust is heated by the interstellar radiation 
field, and a small fraction is heated by stronger radiation 
fields associated with star formation. 

We caution the reader that estimating dust masses 
based on broadband infrared fluxes, as we do here, is 
a poorly constrained technique; there are very few data 
points compared to the number of parameters in the 
model. The DL07 model uses five parameters to charac- 
terize the emission from dust in galaxies: Muust, Umin, 
Umax 1, and a. The dust mass is represented by Muust- 
Umin represents the interstellar radiation field heating 
the diffuse ISM, and Umax represents the upper limit on 
the interstellar radiation field. The starlight heating the 
dust is described using the dimensionless parameter U , 
that by definition is always between Umin and Umax ■ The 
value [/ = 1 is the local interstellar radiation in the Milky 
Way. The parameter 7 represents the fraction of gas 
that is exposed to strong radiation fields with intensities 
in the range Umin <U < Umax- Finally, a characterizes 
the distribution of starlight intensities. In practice we fix 
two of the se parameter s , a an d Umax- We adopt the val- 
ues set bv IDraine k Lil (|2007D of Umax = 10*^ and a = 2. 



Therefore three parameters are free in the model {Muust , 
Umin, and 7). We remind the reader that we constrain 
this model with only four broadband fiuxes, at 8, 24, 70 
and 160 fim. We can also use the returned values of these 
parameters to calculate a temperature for the majority 
of the dust grains {Tjjmin), the fraction of dust luminos- 
ity that originates in photon dominated regions (fpoB.), 
a nd the dust- weighted n iean s tarlight intensity ((C/)). 

iMunoz-Mateos et al.l ()2009D derived empirical fits re- 
lating a grid of DL07 emission model outputs to the 
Spitzer fluxes. In particular, the DL07 parameters 
Moust, 7, (U), and Jpdr can all be derive using mea- 
surements at 8, 24, 70 and 160 /im. We can calculate 
Umin using {U) and 7 according to Eq. (33) in DL07. 
For the I Zw 18 fluxes at 8 and 24 um we used the val- 
ues measured by lEngelbracht et al] (120081 ) of 0.47 and 
6.28 mJy respectively. The agreement between the dust 
mass derived using the empirical flts and DL07 models 
is very good, with a scatter of about 9% and an offset of 
-f 5%. The DL07 dust masses are strongly dependent on 

R70 = {vF^)7o/{i^Fu)i&0: with Moust oc R-TO is 

sensitive to the temperature of the largest grains domi- 
nating the FIR emission, and any new constraint or de- 
tection at 160 //m will strongly affect the resulting dust 
mass. Essentially, the smaller the 160 /xm flux, the hot- 
ter the temperature and thus the less dust is needed to 
produce the observed 70 /im emission. 

The derived dust properties for I Zw 18 are summa- 
rized in Table [ll Median v a lues f or 48 SINGS galax- 
ies analyzed bv lDraine et al.l (|2007[ ) and three starburst 
system out of the same sample are included for com- 
parison. The lower limits obtained for Umin, (U), 7 
and fpDR in I Zw 18 are high co mpared to the mean 
values in the SINGS sample (D raine et al.l 120071 ). The 
high radiation intensity environment of I Zw 18 is com- 
parable to starbursting systems like Mrk 33, T ol 89 
and NGG 3049, as also found bv IWu et al.l (l2Q07l) . We 
flnd that the DL07 model yields a mass upper limit of 
M_D„st < 1.1 X 10^ Mq. Just like the modified blackbody 
case, if we use the MIPS 160 /im upper limit instead of 
the PACS upper limit, we measure a dust mass a factor 
of -2 higher, i.e., Moust < 2.6 x lO'' Mq. 

3.2. Comparison of Dust Masses 

The DL07 dust mass upper limit is a factor ^3.5 larger 
than the dust mass estimated using the modified black- 
body model. DL07 model treats dust emission as an 
ensemble of dust grains at different temperatures that 
includes larger masses of dust at colder temperatures 
than what is predicted by the single temperature fit. 
Therefore, it is not surprising that this model generates 
a higher dust mass than the modified blackbody model. 
However, the fact that these two measurements are not 
extremely different increases our confidence in the dust 
mass limit, which we conservatively take to be that re- 
sulting from the DL07 model. 

It is interesting to compare this result to other mod- 
eling effort s for I Zw 18 and low metallicity galaxies. In 
particular. I Galametz et a l. (2011) determine dust masses 
in a large sample of galaxies with literature data, using 
full spectral energy distribution modelling based on the 
iZubko et ai] ^2004 ) grain model. They find that the in- 
clusion of submm-wave data tends to drive dusty galaxies 
toward lower dust masses, while for low metallicity galax- 
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TABLE 1 

Derived dust properties based on the DL07 model 



Object 


Monst (Mq) 




(u) 


7(%) 


fpDR(%) 


T^Umin (K) 


I Zw 18 


< 1.1 X 10" 


> 8.8 


> 21.1 


> 13.1 


> 48.9 


> 24.4 


Mrk 33 


2.9 X 10" 


4.0 


14.3 


11.8 


47.6 


21.4 


Tol 89 


2.8 X 10^ 


2.0 


3.6 


6.4 


33.2 


19.1 


NGC 3049 


5.5 X 10^ 


3.0 


5.1 


6.1 


32.6 


20.4 


SINGS'* 


1.5 xlO^ 


1.5 


1.7 


0.9 


7.8 


18.2 



Median values for 48 SINGS galaxies (Table 5, [ 



ies the inclusion of submm-wave constraints yields higher 
dust mass p redictions than thos e from far-infrared alone. 
By contrast [Praine et al.l ()2007| ) found their masses to be 
robust to the inclusion of submm-wave data. This is in 
part driven by modeling choices, in particular the inclu- 
sion of a minimum ra diation field or the interpretation o f 
submm-wave excess (|Israel et al.|[2010l: IBot et al.llMol) 
as caused by cold dust. The latter appears to be as- 
sociated to low metallicities, and although re cent studies 
sugge st that it is not caused by cold dust (Ga lliano et al.l 
1201 If) , it remains a large systematic uncertainty in dust 
mass determinations. Also note that f or dwarf gala x- 
ies the SCUBA data included bv Galam etz et al.ll2011lis 
much deeper than the data used by iDraine et al.l 120071 
Observational biases and limitations have a non negligi- 
ble impact on the interpretation of the observed trends 
i n dust-to-gas ratio vers us metallicity. 

iGalametz et alj (1201 ID use a 160 u m flux upper limit 
of 76.8 mJv (jEngelbracht et al.ll2008[ ). For submm -wave 
data they use a 850 /xm flux upper limit of 1.25 mJy 
(IGalhano et al.l 120081 ) . They measure an upper limit to 
the dust mass of 1.1 x 10^ Mq scaled to our adopted dis- 
tance. Our dust mass up per limit is an o r der of magni- 
tude smaller than that of lGalametz et aH (l2011f). Using 
the DL07 model and the data in lOalametz "eTall ([20T1 
we would obtain an upper limit o f 1.0 x 10^ Mp,. This is 
nearly equivalent to the dust mass lGalametz et al.l (|2011l ) 
finds, highlighting the fact that the mass limits are likely 
reasonably robust to the choice of models. Therefore, 
our lower dust mass limit is due to the tighter flux limits 
at 160 /xm. 

3.3. Gas Mass 

The tota l HI mass of I Zw 18 is Mhi = 2.3 x 10^ Mq 
(jvan Zee e t al. 1998). The molecular content of I Zw 18, 
however, remains unknown since no CO emission has 
been detected. Our new upper limit on the CO J = 1 — > 
luminosity of I Zw 18 is Leo < 10^ K km s~i pc^ (4cr), 
which corresponds to Mh2 < 450, 000 Mq for a standard 
conversion factor {aco — 4.5 Mq (K km s~^ pc^)~^). 
Note that our lum inosity is similar to that quoted by 
iLeroy et"al1 ()2007br ) because we adopt here a larger dis- 
tance (for matched smoothing and assumptions we im- 
prove the sensitivity of that study by a factor of 2). Using 
this Milky Way based conversion factor, Mh2 is at most 
0.2% of the total gas mass. 

There is no reason to expect that the Milky Way based 
conversion factor between CO luminosity and H2 mass 
applies to low metallicity galaxies like I Zw 18. In 
the L ocal Group, ctgo is a strong function of metal- 
hcity ()Lerov et all 1201 ID . iGenzel et afl ()2012D derive a 
correlation between oxygen abundance and conversion 



factor, aco- Applying their formula to a galaxy with 
the metallicity of I Zw 18, we find a conversion factor 
aco ~ 477.5 Mq (K km s^^ pc^)"^ This factor is ~ 100 
times larger than aco in a typical spiral galaxy. Using 
this conversion factor we calculate a molecular gas mass 
upper limit of Mm < 4.8 x lO"^ Mq. This Mm is 20% 
of the total gas mass. 

There is some evidence that, at low metallicities, the 
star formation activity may be a better indicator of the 
mole c ular mass than the CO emission (Krumholz et al.l 
[20Tlt iBoIatto et al.l[20nl: ISchruba et al.ll201lD . The Ha 
flux of I Zw 18 suggests a recent star formation rate 
(SFR) of ~ 0.1 Mq yr~^. In large star-forming galax- 
ies, a typic al H2-to-SFR rat io (H2 depletion time) is 

1-2 Gyr (jBigiel et al.ll2011h . The H2 mass correspond- 
ing to this amount of star formation in such a galaxy 
would thus be ~ 10* Mq. Given the level at which star 
formation obviously dominates the morphology and ISM 
conditions in I Zw 18 we strongly suspect that this kind of 
equilibrium assumption very much overestimates the H2 , 
which will have been dissociated or otherwise destroyed 
by the recent burst. Nonetheless, even in this limit the 
II2 only makes up ^ 30% of the integrated gas mass. The 
similarity wi t h the results obtained from applying the 
IGenzel et al.l (|2012[ ) correlation is not surprising, since 
the underlying assumption is the same. Because of its 
uncertainty, we do not include the H2 correction in the 
following calculations. 

3.4. Dust-to-Gas Mass Ratio and Metallicity 

IDraine et all (|2007D finds that the DGR changes sig- 
nificantly depending on whether the dust mass is com- 
pared to the total gas mass or only the gas mass enclosed 
in the aperture where the infrared emission is measured. 
For example, IC 2574 (a dwarf galaxy in the SINGS sam- 
ple), only 19% of the HI gas mass is enclos ed in the area 
where the infrared emission is detected ([Walter et al.l 
120071 ). We take the point of view that these are two 
valid definitions of the DGR: global, or local where the 
dust emission is dete cted. Using the total HI mass from 
Ivan Zee et al.l (|1998[ ). we measure an upper limit for the 
global DGR < 5 x 10~^. Using instead the HI mass 
enclosed in the area were we measure the 160 /im flux 
upper limit (62% of the total HI mass), yields a local 
DGR < 8.1 x 10-^. 

Figure H shows the DGR as a function of oxygen abun- 
dance for I Zw 18 and a subsamplc of SINGS galaxies. 
Open symbols indicate 3(t upper limits. The solid line 
represents a linear scaling of the DGR with metallic- 
ity. This linear relation assumes that the abundances of 
all heavy elements are proportional to the oxygen abun- 
dance and that the same fraction of all heavy elements 



Dust in the Extremely Metal Poor Galaxy I Zwl8 



7 




7 7.5 8 8.5 9 7.5 8 8.5 9 

12 + log(0/H) 12 + log(0/H) 

Fig. 4. — Dust-to-gas mass ratio is shown as a function of oxygen abundance. Open symbols correspond to upper limits (Scr). I Zw 18 
DGR is measured based on the PACS 160 fim flux upper-limit and is shown as an open star. The upper error bar corresponds to the value 
of the DGR if estimated based on the MIPS 160 /im flux (Sir) upper-limit, which may be more appropriate if the source is extended at 160 
^m over scales larger than 12". Circles and triangles correspond to SINGS galaxies with and without SCUBA fluxes, respectively. The 
dust and gas mass values are from Tables 4 and 5 in Draine et al. ( 2007). The values in the left panel represent the global DGR when 
estimated using the total gas mass of the galaxy. The right panel shows the DGR values when they are estimated "locally" , using only the 
gas mass in the region where the IR em ission is detected. The solid line shows a linear scaling between DGR and metallicity normalized 
to match the Milky Way values (eq. 13. IDraine et al.ll2007l '). 



are i n solid form as in the Milky Way (jPraine et alj 
I2007D . The I Zw 18 DGR upper limit is primarily 
driven by the upper limit in the dust mass, while in 
the SINGS galaxies the upper limits are due to lower 
limits in the gas mass (due to the non inclusion of 
H2). As we discussed in Section 3.1, we obtain differ- 
ent dust masses for I Zw 18 depending on the assump- 
tion we make about the distribution of the 160 /im emis- 
sion (point-like with PACS and extended with MIPS). 
The open star shows the DGR of I Zw 18 when we as- 
sume point-like emission, while the upper limit of the 
bar shows the DGR when we assume extended emis- 
For the SINGS g alaxies , the dust and gas masses 

and the me tallicities are 
Note that iDraine et al.l 



sion. 



are f rom IDraine et aTT (12007 
from Moustakas et al.l (I201C 



()2007| ) computes molecular gas masses assuming a fixed 
Xno factor of 4 x 10^° cm~^ (K km s~^)~^ following 
(jBlitz et al.ll2007| ). We show SINGS galaxies with and 
without measur ed SCUBA fluxes as triangles and circles 
respectively, as IDraine et al.l ()2007D find that the dust 
mass estimates with and without SCUBA data can dif- 
fer by a factor of ^2. 

The left panel of Figure H] shows the global DGR, es- 
timated using the total gas masses. The right panel 
shows how the DGR changes when estimating it locally 
in low metallicity systems, including I Zw 18, by using 
the gas mass enclosed in the region where the infrared 
emission is detected. For SINGS galaxies with metal- 
hcities 12 -I- log(0/H) > 8.1, the total DGR seems to 
agree within a factor of ^2 with a linear relationship be- 
tween DGR and metallicity. Low metallicity galaxies do 
not seem to follow the same linear correlation that in- 
cludes the Milky Way DGR. The I Zw 18 global DGR 
falls below the linear scaling by a factor of ^ 8. The 



right panel of Figure |4] shows the local DGR. The DGRs 
of the SINGS low metallicity systems scale up and appear 
consistent with the linear relationship within a factor of 
'^2, although most of the low metallicity points are only 
upper limits. For I Zw 18, however, the local DGR falls 
below the linear scaling by a factor of ^ 5. Therefore, 
our dust mass limits for I Zw 18 suggest a breakdown 
of the linear relationship between DGR and metallicity 
at very low metallicities. Note further that, at least in 
terms of the global DGR, I Zw 18 seems to continue the 
trend found for other low metallicity galaxies. 

We show in Fig. |4] that only one of the seven SINGS 
galaxies with 12 -I- log(0/H) < 8.1 has a DGR that 
is not an upper limit. It may be possible that the lo- 
cal DGR of this system is higher than other low metal- 
licity galaxies and the trend is really steeper than lin- 
ear, as our result for I Zw 18 and other studies suggest 
(jLisenfeld fc Ferraral 119981: iMufioz-Mateos et al.l [2009il . 
Clearly more work is needed to determine, robustly, 
wether low metallicity galaxies do or do not follow the 
linear scaling shown in Fig. (4] 

Note that the abundance of oxygen may not be the 
correct abundance to refer to. Indeed, the abundances 
of refractory elements that constitute the bulk of the dust 
such as carbon or silicon are like l y mor e relevant to estab- 
lishing the DGR. iGarnett et al.l ()1999l ) , for example, find 
a trend of increasing C/0 with O/H for a sample of irreg- 
ular and spiral galaxies observed with HST. This could 
suggest that the nonlinear trend of DGR with metallicity 
is really an artifact of using O/H as a proxy for metal- 
licity, and the relation coul d become more nearly linear 
when plotted against C/H. IGarnett et al.l (|199 9^ find a 
gas-phase abundance of C in I Zw 18 that is significantly 
higher than that predicted by the extrapolation of the 
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observed C/0 vs. 0/H trend in low metallicity irregular 
galaxies. In fact, C/0 in I Zw 18 is only 0.3 dex lower 
than Solar. This is barely enough to reconcile our lim- 
its on the local DGR with a linear trend with C/H, and 
probably not enough to explain our low global DGR, but 
it certainly goes in the right direction. 

4. CONCLUSIONS 

In this work we study I Zw 18 using data from Spitzer, 
Herschel Space Telescope and IRAM Plateau de Bure In- 
terferometer. We reduce the flux upper limit at 160 fim 
by a factor of ^3 and the CO J = 1 ^ flux up- 
per limit by a factor of '^2 compared to previous mea- 
surements. Combinin g these observation s with the dust 
emission model from iDraineini (|2Q07l ). we constrain 
the dust mass to be Mdust < 1.1 x 10^ Mq. We note 
that any dust mass measurement relies on assumptions 
about the mass emissivity of dust grains in the inter- 
stellar medium, with the important associated system- 
atic uncertainties. We find a global dust-to-gas mass 
ratio of Mdust/^gas < 5.0 x 10^*^, while the ratio mea- 
sured in regions where the 70 /im emission peaks is 
Mdust/Mgas < 8.1 X 10-5. 

These measurements are suggestive that low metal- 
licity galaxies do not follow the same linear relation- 
ship between metallicity and DGR as typical local spi- 
rals. At face value our DGR upper limit is inconsistent 
with the hypothesis that the fraction of heavy elements 
incorporated into dust is the same in high metallicity 
galaxies (such as the Milky Way) and in extremely low 



metallicity galaxies (such as I Zw 18). There are other 
scenarios, however, that can produce a break or non- 
linear power-law relationship between DGR and metal- 
licity. For instance, models that include more detailed 
physical processes such as the production and destruc- 
tion of dust by supernovae, removal of dust through 
outflows from galax ies, and dust production in the en- 
velopes of stars (e.g.lLisenfeld &: Ferraralll998tlEdmundsl 
l200lt IHirashita et al.l I2002D may yield non-linear rela- 
tions. Much more work is needed with sensitive maps 
of low metallicity galaxies, like I Zw 18, to better under- 
stand the realtionship between DGR and metallicity. 
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